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Double quantum dot as a spin rotator 
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It is shown that the low-energy spin states of double quantum dots (DQD) with even electron 
occupation number N possess the symmetry 50(4) similar to that of a rigid rotator familiar in 
quantum mechanics (rotational spectra of H2 molecule, electron in Coulomb field, etc). The "hidden 
symmetry" of the rotator manifests itself in the tunneling properties of the DQD. In particular, the 
Kondo resonance may arise under asymmetric gate voltage in spite of the even electron occupation 
of the DQD. Various symmetry properties of spin rotator in the context of the Kondo effect are 
discussed and experimental realization of this unusual scenario is proposed. 
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I. INTRODUCTION 



During recent years, the physics of single electron tun- 
neling through a quantum dot (QD) under the condi- 
tions of strong Coulomb blockade has been at the fo- 
cus of intense investigation 111 . The number of electrons 
N in a dot can be regulated by a suitable gate voltage 
V g applied to an electrode coupled capacitively to the 
dot. The Coulomb blockade suppresses the tunneling 
through the dot unless the resonance between its ener- 
gies filled by N and N + 1 electrons occurs at certain 
values of V g , when it compensates the charging energy, 
i.e. £(N + l,V g ) ~ £(N,V g ). The differential conduc- 
tance dl /dVsd of a QD forms a diamond-like patterns in 
the plan (V s d,V g ) where the non-conducting "windows" 
are separated by a network of Coulomb resonance lines 
(here V s d is the source-drain voltage) . 

Accurate low-temperature experiments demonstrated 
the existence of Kondo resonances in the windows corre- 
sponding to odd occupation of the dot |g{ (O-diamonds). 
These resonances are seen as zero-bias anomalies (ZBA), 
i.e. as bridges of finite conductance connecting two op- 
posite vertices of O-diamond-shape window at V s d — ► 0. 
Besides, it was predicted theoretically || and observed 
experimentally Q that Kondo resonances can appear also 
in the even occupation windows (E-diamonds) at strong 
enough magnetic fields. This unconventional magnetic 
field induced Kondo effect arises because the spectrum of 
the dot possesses a low-lying triplet excitation when the 
electron at the highest occupied level is excited with spin 
flip. The Zeeman energy compensates the energy spac- 
ing between the two adjacent levels, and the lowest spin 
excitation possesses an effective spin 1/2, thus inducing 
a Kondo-like ZBA in the differential conductance. 

Similar effect is possible in vertical quantum dots for 
which the singlet and triplet states may be close in energy 
both at even and odd occupation. The influence of an ex- 
ternal magnetic field on the orbital part of the wave func- 
tions of electrons in vertical quantum dots is, in general, 



more pronounced than the Zeeman effect. Hence, singlet- 
triplet level crossing are induced by this field, causing the 
emergence of Kondo scattering at even filling or its en- 
hancement at odd filling . The theory of Kondo tunnel- 
ing through vertical quantum dots in external magnetic 
field is developed in Refs 



In the present paper we explore yet another device 
which manifests the Kondo effect in QD with even elec- 
tron number N, namely, a QD with two wells which is 
refereed to as double quantum dot (DQD). A systematic 
treatment of the physics of DQD with even N coupled 
to metallic leads is presented below. Special attention is 
given to the symmetry properties of DQD and its repre- 
sentation as a quantum spin rotator. It is well known 
that the tunneling Hamiltonian for a QD can be mapped 
on the Kondo Hamiltonian in the O-diamond window of 
the QD. In the E-diamond window, the same procedure 
of eliminating the charged virtual states results in the 
four-state Hamiltonian of doubly occupied dot where the 
singlet 5 = and triplet 5 = 1 levels are intermixed 
by second-order tunneling. As is shown in Ref. |10[, this 
effective Hamiltonian possesses the dynamical symmetry 
50(4) of a spin rotator. As a Kondo scatterer, spin ro- 
tator possesses new properties in comparison with local- 
ized spins obeying SU(2) symmetry. The magnetic field 
induced Kondo effect mentioned above is a manifestation 
of a " hidden symmetry" which is a footprint of the 50(4) 
group. 

In Ref. [|l(J a special case was considered, namely, an 
asymmetric DQD formed by two dots of different radii 
in a parallel geometry coupled by tunneling interaction 
with even occupation N = V\ + v r (I, r stands for left and 
right respectively). Moreover, it has been assumed that 
the strong Coulomb blockade exists in one dot whereas 
tunneling contact with the metallic leads exists in the 
other dot. Here we will address more general situations 
and compare several representations in terms of effec- 
tive spin Hamiltonians. It will be shown that unusual 
ZBA can arise in generic DQD structures. In particular, 



the Kondo effect induced by quantum dots with SO (4) 
spin rotational properties exists in asymmetric DQD also 
when both the I and r dots are coupled with the leads 
and the Coulomb blockade is strong enough in both of 
them. The main precondition for the emergence of the 
Kondo effect in this case is the sizable difference in ion- 
ization energies of the two dots. This quantity can, in 
fact, be tuned by an application of a suitable gate volt- 
age to one of them. The same effect can be achieved also 
in a symmetric DQD with even occupation in a parallel 
geometry provided the axial symmetry of the system is 
broken by the difference in gate voltages applied to the 
right and left dot (Vp 1 , respectively). 

DQD oriented parallel to the lead surfaces were fab- 
ricated several years ago OJl2]. Two main resonance 
effects were noticed in such electric circuit. First, one 
of the dots (say, right) can be used as an electrome- 
ter p2| . Scanning V£ at fixed Vg, Coulomb oscillations 
can be induced both in the right and left dot because 
the interdot capacitive coupling changes the positions 
of Coulomb resonance in both of them. As a result, 
the step- wise structure of the conductance acquires more 
complicated form. The Coulomb blockade windows be- 
tween the resonances in the Coulomb energy of the dot 
£v r ,vi(V£,Vg) form an "egg-carton" pattern Q] where 
the vertices connect the windows with charge configu- 
rations (y r ,i/{),(v r ,vi — 1), [v r + l,Pi — I). The lines 
£i/j,i/ r ~ £u t +iv r , are the regions where the Coulomb res- 
onance induced by VT allows tunneling through the left 
dot. Second, the resonance £ Vr>ul ~ £„,.+!,„, _i allows co- 
tunneling through the right and left dots, which is the 
precondition for the Kondo effect due to appearance of 
pseudospin-like configuration of the DQD jD|. Then, 
manipulating with VI, one can induce the third transi- 
tion i/i — l,v r -\-l — ► Ui, v r thus closing the loop and orga- 
nizing the "electron pump" which transfer single electron 
from one dot to another (see |1J] and references therein). 

The picture becomes even richer if the tunneling be- 
tween the right and left wells of the DQD is taken into 
account. Then the dot can be treated as an artificial 
molecule where the interdot tunneling results in for- 
mation of complicated manifold of bonding and anti- 
bonding states Jl5[ which modifies its charge degrees 
of freedom. Besides, it induces an indirect exchange, 
thus modifying the Kondo resonances when the dots are 
placed in series m&. ft will be shown below that the 
interdot tunneling in parallel geometry results in the ap- 
pearance of Kondo precursor of Coulomb resonance along 
the lines £ Vr ., Vl ~ £v r +\,vi-\ provided there exists di- 
rect tunneling coupling V between the left and right 
dot. We consider the simplest case of v\ = v r = 1 in 
a neutral ground state of DQD. ft will be shown that 
unconventional Kondo resonance occurs under condition 
V/[£i t 2 — £1,1] <C 1. Moreover, this kind of Kondo res- 
onance can appear also in the middle of the Coulomb 
window for the right dot, provided the capacitance of 



the left dot exceeds essentially that of the right dot jf?] . 
In both cases, the DQD possesses the symmetry of a spin 
rotator. 

In section 2 the various setups of DQD are introduced 
and the Hamiltonian describing the DQD is written down 
within the framework of a generalized Anderson model. 
The phase diagram of charging states in the left and right 
gate voltages plan is schematically drown, and the regions 
of Kondo resonance are indicated. In the first part of 
Section 3 the spectrum of the isolated dot with even oc- 
cupation is discussed. The second part is devoted to the 
derivation and solution of renormalization group (RG) 
equations for the DQD. The central result of this sub 
section is a demonstration of possible singlet triplet level 
crossing due to tunneling. When the renormalized en- 
ergies are below the reduced band edge, renormalization 
stops, and charge fluctuations are suppressed. This is 
the Schrieffer- Wolff regime, and a derivation of an effec- 
tive spin Hamiltonian is executed in section 4. In the first 
part, the spin Hamiltonian is given in terms of two vector 
operators and is shown to have the SO(A) symmetry of a 
spin rotator. It is followed by a short subsection in which 
the renormalization group flow of coupling constants is 
explained and the Kondo temperature is derived. Then, 
in the third subsection, a two spin representation is sug- 
gested, in which the occurrence of two spin 1/2 operators 
just reflects the fact that the algebra 04 is a direct sum 
of two 03 algebras. In the fourth subsection, the possi- 
bility of arriving at the Kondo effect in finite magnetic 
field is discussed, leading to a third representation of the 
spin Hamiltonian. The question of whether a DQD with 
two electrons can be regarded as a real two site Kondo 
system (even if the DQD is highly asymmetric) is dis- 
cussed in section 5. In particular, a stringent comparison 
is made with the two spin representation mentioned in 
section 4. The paper is concluded in section 6. Some 
technical details of various calculations are relegated to 
the Appendix. 



II. MODELS OF DOUBLE QUANTUM DOTS 
WITH SINGLET GROUND STATE 

Two models considered in this work are sketched in 
Fig.l. We will refer to a system (a) with zero gate volt- 
ages as a "symmetric" DQD. The same system with Unite 
but unequal gate voltages V g ' r will be called a "biased" 
DQD and the pair of dots with different radii shown in 
Fig. lb will be referred to as "asymmetric" DQD. 




(a) 



left 



right 




FIG. 1. Double quantum dots in parallel geometry. Left (I) 
and right (r) dots are coupled by tunneling V to each other 
and by tunneling Wi iT to the source (S) and drain (D) elec- 
trodes. V g ' r are the gate voltages, (a) symmetric dot, (b) 
asymmetric dot. 

In all cases the DQD is described by a generalized An- 
derson tunneling Hamiltonian which takes into account 
the internal structure of the DQD, 

H = H b + H t +H d + H g . (1) 

The first term, Hj, is related to the lead electrons, 

ka,a 

where a — s,d denotes electrons from the source and 
drain electrodes pa] and a = ± is the spin index. The 
second term, H t is the tunneling Hamiltonian, 



H t =J2J2 (wwcLdto + h.c.) 



(3) 



i—lr ka 



Here c ka = 2~ l / 2 { 

w ka ,i/(wi s4 + W 2 d - 



Cka.s + Ck<r,d), and W k i = 
W kd,i) 1/2 - Thc third term > H dl describes 
the isolated DQD. In the present context, the quantum 
dot is a " molecule" , containing N — viq + v r $ electrons 
in a neutral ground state. The capacitive interaction be- 
tween the two wells of the DQD is assumed to be strong 
enough to suppress the fluctuations of electron tunneling 
induced occupation in the windows between the Coulomb 
resonances of tunneling amplitude. We consider DQD 



with even TV, so that, generically, the ground state of an 
isolated DQD is a spin singlet. The isolated dot is then 
described by the Hamiltonian, 



H d = ^2 ^2 e i n ia + V ^ d L d J° + H corr, (4) 

i=l,r (7 i^j 



in which V is the inter-well constant tunneling amplitude. 
The capacitive interaction within the DQD is described 
by the term 



H,.. 



5?' 



ini(n t - 1) + Qi r 5nidn r . (5) 



Here n, — Y2 a d i(J di ai and 5rii = i%i — ViQ is the deviation 
of electron distribution from the neutral charge config- 
uration v M for a given DQD. Moreover, Qi — e 2 jld is 
the charging energy of the dot i whose capacitance is Ci , 
and Qi r is the capacitive coupling between the left and 
right dots. The simplest configuration which contains 
in a nutshell all the complicated physics of many-body 
interactions arising in a course of tunneling is N = 2, 
f/o = v t q = 1. This case, for which Srii = rii — 1., will 
be given a special attention below. Finally, the term H g 
represents the gate voltage energy. We consider symmet- 
ric and asymmetric DQDs formed by wells of equal and 
different radii respectively (Fig. 1). Hence, generically, 
the gate potential H g is asymmetric, 



tf 9 = E^ v l^ v l 



(6) 



It is, in fact, useful to include the gate potential (^) in the 
position of the one-electron energy levels, £j = e^ + V g . 
Then, by tuning the gate voltage one can change the en- 
ergy difference A = Si—e r or, in other words, redistribute 
the electron density between the left and right wells of 
the DQD. 

It is assumed that in equilibrium and at zero gate volt- 
ages, each dot is filled by one electron and the Fermi level 
of the leads is in the middle of the Coulomb blockade win- 
dow. The energy levels of a symmetric DQD with uncou- 
pled dots (Qi = Q r = Q, V = 0)) are shown in the upper 
panel (a) of Fig. 2. These levels may be shifted relative 
to each other and to the Fermi level £f, and each level 
crossing Ei — Ef corresponds to recharging of the dot i. If 
electron exchange between the right dot and the leads is 
blocked M], the charge transfer resonance between the 
states {1, 1} and {0, 2} occurs when Ei = £ r + Q (see also 
fll9|). In the general case (Fig. la), additional electron 
appears in the dot i when the levels £i + Q and Ef cross. 
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FIG. 2. Energy level scheme for symmetric (a), biased (b) 
and asymmetric (c) DQD. Upper panel: filled and empty 
one-electron levels. Dashed arrows indicate charge-transfer 
excitons. Lower panel: two-electron states of isolated and 
coupled left and right dots. 

In the absence of interdot tunneling, V = 0, one can 
easily obtain the effective spin Hamiltonian for the DQD 
with TV = 2 in the ground state far deep in the Coulomb 
blockade windows. This is the two-site Kondo Hamil- 
tonian in the window {1,1} and the single site Kondo 
Hamiltonians in the windows {2,1} and {1,2}. In the 
latter case of charged DQD occupied by odd number 
of electrons, tunneling through the left (right) dot is 
blocked, but a Kondo-type resonance compensates for 
the Coulomb blockade and opens a tunneling channel 
through the right (left) dot. In the former case of neu- 
tral DQD with even occupation the possibility of Kondo 
tunneling is determined by the relative strength of the on- 
site indirect exchange Ji between the spins Si of singly 
occupied dots and conduction electrons in the reservoir 
on the one hand, and the sign and magnitude of the in- 
tersite RKKY exchange Ji r on the other hand p(J . Both 
these parameters are predetermined by the tunnel cou- 
pling constants Wki with the band electrons in the reser- 
voir, but one can modify them by varying the gate volt- 
ages and interdot distance. 

The interdot coupling significantly modifies this pic- 
ture. It favors the singlet spin state in the middle of 
the Coulomb blockade window {1,1} thus eliminating 
the Kondo tunneling at zero gate voltages. At finite 
\Vl — Vg\ the values of v\ deviate from the integer val- 
ues near the boundaries between the different charge sec- 
tors. Increasing negative gate voltages V* or VJ, one 
can bias the charge distribution in favor of left or right 
dot, respectively, without changing the total number of 
electrons. As a result, with increasing \Vg — V g r \ one 
reaches the region of states with small charge transfer 
gap Ai = Q + e r — Si <C Q. The energy levels of such 
"biased" DQD are shown in the upper panel of Fig. 2b. 
These states occupy the upper right corner of the window 
{1, 1} hatched in fig. 3a. Here the zero energy configura- 



tion illustrated by Fig. 2a corresponds to the coordinate 
origin. The virtual charge transfer excitations (dashed 
arrow in the upper panel of Fig. 2b) significantly influ- 
ence the tunneling through the DQD. It will be shown 
below that a novel type of Kondo resonance arises in 
this area of the sector {1, 1}. The "biased" DQD in this 
sector behaves like a spin singlet at high temperatures 
and excitation energies, and demonstrates the properties 
of spin one triplet partially screened by the Kondo tun- 
neling at low energies and temperatures T < Tk- The 
Kondo temperature Tk is a function of V* . The Kondo 
"isotherm" TK(Vg,VI) is presented by the dashed line in 
Fig. 3a. 
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FIG. 3. Coulomb windows for different charge states 
{yi, Vr\ of symmetric (a,b) and asymmetric (c) DQD. Hatched 
regions indicate the domains where the Kondo effect exists. 

Similar effect exists for the asymmetric DQD (Fig. lc) 
where the two coupled dots have different radii r\ ^> r> 
and hence different blockade energies, Qi -C Q r . The 
energy levels of an isolated DQD are shown in the up- 
per panel of Fig. 2c, and the corresponding recharging 
map is presented in Fig. 3b. Here the hatched area also 
marks the region of the map where Kondo effect arises 
in spite of the even number of electrons in the dot. The 
Kondo isotherms in this case are parallel to the Coulomb 
resonance line. 

Conventional approach for the description of the 
Kondo effect in a two-site quantum dot starts with the 
two- center Hamiltonian 



Hdo = Hi + H r , 



(7) 



and treats H t in terms of a two channel tunneling oper- 
ator, 



Ht — H t i + H t 



(8) 



The interdot interaction Hi r is considered as a coupling 
between two resonant Anderson centers. If the left and 
right dot each contains an odd number of electrons (as in 
our simplified model with v^ r = 1), the Kondo tunneling 
is possible through each dot separately. The exchange 
part of interdot coupling maps our Hamiltonian onto the 
two-site Kondo model. This coupling can be both of fer- 
romagnetic and antiferromagnetic type. In the latter case 
the interplay between H t and Hi r results in suppression 
of Kondo tunneling through the left and/or right well 
of the DQD. The phase diagram of the two-site Kondo 



model is discussed in numerous papers |20|. 

In the model discussed here, the interdot interaction 
is represented by the term Hi r = V^,- 9 t,-d] (7 .dj <T till we 
remain in the charge sectors {1, 1}, {1, 2}, {2, 1} of Fig. 
3. It is obvious that this coupling suppresses Kondo tun- 
neling through the symmetric DQD at zero gate volt- 
ages (point "0" in Fig. 3a) because the effective in- 
direct exchange interaction which arises due to virtual 
excitations of charged states {0, 2}, {2, 0} is of antiferro- 
magnetic sign, Ji r — 2V 2 /Q like in the Heftier-London 
limit for a hydrogen molecule or in the half-filled Hub- 
bard model. As a result, the ground state of a DQD is a 
spin singlet, and the gap S = Et—Es = Jw which divides 
the triplet excitation from the singlet spin ground state 
prevents the formation of a Kondo resonance. It will 
be demonstrated below that this is not so in the case of 
strongly asymmetric DQD (hatched regions in Fig. 3a, b), 
where the crossover to a triplet state is induced by the 
tunneling H t . 

To describe this crossover, it is more convenient first 
to diagonalize the dot Hamiltonian H c i, i.e. to express it 
in the form 



H d = Y,E N A\NA)(NA\, 



(9) 



JV,A 



recalling that N is the number of electrons in a given 
charge state of the DQD whereas A stands for a set of 
quantum numbers which characterize the many-electron 
configuration d" 1 d v T r in the presence of interdot coupling. 
In order to get compact form for some equations, we 
introduce Hubbard projection and configuration change 
operators 

X NA > N ' A ' = \NA)(N'A'\. (10) 

The diagonal terms X NA ' NA are conventional projection 
operators, while the off-diagonal operators change elec- 
tron configuration of the dot. The tunneling term H t (g) 
can now be rewritten in the form, 



Hi 



EEE 

JV,AiV',A' ka 



(lTrNh,N'A' Y NA,N'A'. 



h.c. 



(11) 



The matrix elements W. ' are nonzero for states 

in adjacent charge sectors of the eigen space of Hd, 



so that N — N' + 1. In this approach, the DQD is 
treated as a "resonance impurity" in the framework of 
the conventional Anderson model, and its specific fea- 
tures are manifest in a characteristic energy spectrum 
En A which includes contributions due to interdot tun- 
neling and Coulomb blockade. 



III. DOUBLY OCCUPIED DQD AS AN 
ANDERSON IMPURITY 



A. Energy levels and wave functions 

Let us now employ the above approach to the doubly 
occupied DQD with N = 2 in a charge sector {1,1} of 
the Coulomb blockade diagram (Fig. 3) . The dot Hamil- 
tonian (§]) , (|Jj) , (|J) can be exactly diagonalized by using 
the basis of two-electron wave functions 

l 5 > = ;4l>4A|o> ; (i 2 ) 

i f °> = 4|E44 s |o), io = 44Jo> ; 

lea*) = «!<)>, K,) = 4 T 4J0). 

Generically, the spectrum of a neutral DQD consists of 
a singlet ground state 15), a low-energy spin-one triplet 
exciton |T/i) and two high-energy charge-transfer sin- 
glet excitons \Exi) and \Ex r ). The corresponding two- 
electron wave functions are the following combinations: 

\S) = a ss \s) +a s i\exi) + a sr \ex r ) ; (13) 

|T0) = |*o>, \T±) = \t±) ; 
\Exi) = a u \exi) + ai s \s), \Ex r ) = a rr \exi) + a rs \s). 

In the special case of symmetric DQD (e; — e r , Qi = 
Q r ), the axial symmetry allows one to introduce even (e) 
and odd (o) excitonic states 

\ex e>0 ) = -i=(\exi) ± \ex r )). 

The interdot tunneling leaves intact odd singlet state 
\ex ) as well as odd triplet states \T/j,). As a result, one 
has instead of (JO), 

\S) = a ss \s) + a se \ex e ), |7» = |i M ), 
\Ex e ) = a ee \ex e ) + a es \s), \Ex ) = \ex ), (14) 

(see Fig. 2a). We are mainly interested in the limit- 
ing cases of strongly biased symmetric DQD where the 
interdot tunneling results in sizable charge transfer be- 
tween the left and right dot with charge transfer energy 
Ai = e r + Q — ei (Fig. 2b) and asymmetric dot with 
charge transfer energy A2 = ei + Qi — e r (Fig. 2c) . The 



virtual charge transfer transitions which contribute to the 
lowest part of the energy spectrum are marked by the 
dashed arrows. Charge fluctuations are negligible when 

= V/Q <C 1, X = V/Ai <C 1, 02 = V/A 2 < 1, 

(15) 

in cases (a), (b), (c) respectively M]. The expansion 
coefficients a^ in this limit are calculated for all three 



cases in the Appendix (eqs. Al,A2,A3, see also com- 
ment and references [|22|). In a symmetric configuration 
(a) the two-electron levels which correspond to the bare 
states of a symmetric DQD form a low-energy quartet 
w a> t — 2e, uiex e ,o = 2e + Q. The odd states remain un- 
renormalized as a result of interdot tunneling, whereas 
the even states undergo a level repulsion. In the limit of 
small = V/Q < 1 

E s = 2e- 2V0, E T = 2s, (16) 

E = 2e + Q, E e = 2e + Q + 2V0, 

(see lower panel of Fig. 2a). In the case (b) of biased 
DQD, the two-electron bare energy levels are arranged 
as is shown in the lower panel of Fig. 2b: 



u> a ,t = £i 



= ev 



v,l = £l + Q 



The parity is broken in this case, and the singlet state 
\s) is now hybridized with both excitonic states. In the 
limit (051) one has, 



E s -- 


= ei+e r -20 1 V, E T =ei+e r , 


Er-- 


= 2e r + Q + 20 X V, E L = 2e t + Q 



(17) 



20[V 



The role of the "left" exciton \El) in the low-energy pro- 
cesses is negligibly small. In the same spirit, all terms 
~ 0[ — V/(ei — e r + Q) will be neglected in pertinent 
calculations below. 

The general scheme of energy eigenvalues in case (c) 
of an asymmetric DQD is similar to that of case (b), but 
here, the first singlet excitation state is a charge-transfer 
exciton in the left dot. So the bare two-electron spectrum 
is Lu s . t — ei + s r , u> ex j = 2ei + Qi, uj ex ^ r = 2e r + Q r , and 
the hybridized states are approximately given by, 



E s -~ 


= ei+e r - 20 2 V, E T = sl+ s r , 


E L - 


= 2e h + Qi + 2/3 2 y, E R = 2e r + Q r 



(18) 



W 2 v, 



(see lower panel of Fig. 2c). In this case we neglect 
the contribution of "right" exciton \Er) and all terms 
~ /% = V/(e r - £, + QA 

It is seen from Eqs ([16]), (p7|), ( ^8[ ) for the energy spec- 
trum of an isolated DQD that the low-energy excitations 
with energy 8 a = Et — Eg are dominantly of spin char- 
acter, whereas the charge excitations E e<0 in case (a) and 
El.r in cases (b,c) are separated from the ground state 



by the gaps 8 c h ^> 5 S in all three cases under consid- 
eration (lower panels of Fig. 2a-c). This same kind of 
"spin-charge separation" persists when the DQD is hy- 
bridized (via H t ) with itinerant electrons in the metallic 
reservoirs, on which we now focus our attention. 



B. Renormalization of energy levels 

The spectrum of electrons in the reservoirs is contin- 
uous and form a band with bandwidth 2Dq. In accor- 
dance with the renormalization group (RG) procedure 
widely used in the conventional Anderson model, the low 
energy physics can be exposed by integrating out the 
high-energy charge excitations in a framework of poor 
man's scaling technique |23| . This procedure implies 
renormalization of the energy levels and coupling con- 
stants of the Hamiltonian (|l|) by mapping the initial en- 
ergy spectrum — Dq < e < Dq onto a reduced energy 
band -D + \SD\ < s < D - \SD\. 

The mapping procedure results in the following equa- 
tions for the singlet and triplet renormalized energies of 
the DQD: 



E A ^E 



,(0) 



(0) 



EE 



\w. 



AA|2 
q<7 I 



E A 



Ex 



(19) 



where E A is the energy before renormalization, q — 
q u , qb are electron momenta such that e q belong to the 
layers \SD\ near the top or the bottom of the conduction 
band respectively. They appear as intermediate virtual 
states in the processes of positive and negative ioniza- 
tion of the DQD. The index (A = S,Tfj.) in these equa- 
tions is reserved for the neutral two-electron states ( |l3| ) of 
the DQD, whereas the positively and negatively charged 
states with one and three electrons are designated by the 
index A. The wave functions and energy levels of these 
states as well as the matrix elements W^ * are calcu- 
lated in the Appendix. Figure 4 illustrates the processes 
involved in the level renormalization in all three cases 
under consideration. Note that these RG equations are 
uncoupled in this order. In accordance with the poor 
man's scaling approach |23| only the virtual transitions 
with energy ~ D are relevant, and the estimate of the 
sum in the r.h.s. of eq. ( |19| ) gives 



E A = E^- T -^l. 

A A D 



(20) 



where T A = 7rpo |T / ^ A | 2 , po is the density of electron states 
in the reservoir, which is taken to be constant, and W A 
are effective tunneling matrix elements calculated in the 
Appendix. 

The crucial difference between the symmetric config- 
uration (a) and the asymmetric configurations (b,c) is 
that the tunneling amplitudes of the processes involved 
in renormalization (3.1a) are different for singlet and 



triplet states of the DQD. In the symmetric case (a), 
the left and right dot states are involved in renormaliza- 
tion of the two-electron states on an equal footing. The 
relevant processes are \S) — > \eq u , le), 15) — ► \oq u , lo), 
\T) —> \eq u , le), \S — > eq u ,le). The one-electron tunnel- 
ing transitions that give dominant contribution to these 
processes are shown by the dashed arrows in Fig. 4a. 
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FIG. 4. The particle states, which are removed from 
half-filled conduction band on reducing the bandwidth by 
\8D\. The one-electron levels renormalized as a result of this 
process are shown by bold lines, (a) symmetric DQD, (b) 
biased DQD, (c) asymmetric DQD. 

As a result, the tunneling rate in this case is 

r s = r T ^irp Q (\w l \ 2 + \w r \ 2 ). (21) 

In the asymmetric configurations (b) and (c) the even- 
odd symmetry is broken, and the Coulomb blockade 
in one center controls the tunneling through the other 
one |[Ll| , [l2| . Processes with energy ~ D involve only 
the electrons from the left dot. In case (b) the rele- 
vant processes are |5) — ► \q u a,lba), \T0) — ► \q u a,lba), 
\T±) —* |g u ±, lfc±), and the tunneling transitions, which 
give the dominant contribution to these processes are 
shown by the dashed arrow in Fig. 2b. The same kind 
of asymmetry takes place in case (c) (dashed arrows in 
Fig. 4c). As a result, one has, instead of (pi|), 



r T « Tr Po \Wi\ 2 r s = a 2 as r T . 



(22) 



Here the coefficient a ss < 1 is a measure of charge trans- 
fer from the left dot to the right dot due to admixture of 
singlet excitonic states to the ground state singlet [see 
Eqs. (A. 11), (A. 12) in the Appendix]. Iterating the 
renormalization procedure (vw, one comes to the cou- 
ple of differential scaling equations 



dE A _ r A 
dlriD ir ' 



T,S, 



(23) 



which describe the evolution of the two-electron energy 
states with reducing the energy scale of the band con- 
tinuum. These equations describe not only the renor- 
malization of the low-energy two-electron spin states 
but also the change of the one-electron transition ener- 
gies E A — E\, because the one-electron states Ea=i& = 
£ r — 0(/3) are deep under the Fermi level, and the reduc- 
tion of the energy scale does not influence them. 

Scaling equations of the type (p3|) were analyzed in Ref. 
JlO| for a specific case of " Fuldc molecule" or double shell 
quantum dot (DSD), where the electrons in one shell are 
subject to strong correlation effect (Coulomb blockade) 
whereas the loosely bound electrons in the second shell 
are responsible for tunneling, and tunneling to the leads 
is allowed only for the second shell. The model (c) is a 
natural extension of DSD because in the lowest approxi- 
mation in the interdot interaction the tunneling through 
the right dot gives no contribution to level renormaliza- 
tion. In case (b), both left and right dot contribute to 
the renormalization procedure, but the crucial property 
of scaling equations, T$ < Tt (see eq. g2j) is shared by 
both configurations. 

The scaling invariants for equations (E3h are 



7T V 1 A 



(24) 



Here the scaling constants have to be chosen to satisfy 
the boundary condition E a (Dq) = E A . Due to rela- 
tion ( p2| ) the energy Et(D) decreases with D faster than 
Es(D), so that the two scaling trajectories E A cross at 
a a certain bandwidth D = D r estimated as 



r T 



zl± ln Rl =E ^-Ef^5 Q . 



(25) 



According to calculations performed in Ref. pi, this 
level crossing can occur either before or after the 
crossover to the Schrieffer- Wolff regime when the one- 
electron energies E A (D) — En, exceed the half-width of 
the reduced continuous spectrum D ~ \E A (D) — Eu\. In 
both cases, the charge degrees of freedom arc quenched 
for excitation energies within the interval —D<e<D, 
and Haldane's renormalization procedure should be re- 
placed by the Anderson poor man's scaling p3j. 



IV. SPIN HAMILTONIAN FOR DQD 

A. The quantum rotator representation 

The Schrieffer- Wolff transformation [E5[ for the con- 
figuration of two electron states of a DQD projects out 
those states of the dot having one or three electrons and 



maps the Hamiltonian H onto an effective spin Hamil- 
tonian H acting in a subspace of two-electron configura- 
tions (A\...\A'\. 



~ ( — 1 \ m 

H = e s He- s = H + V ^j-[S, [S... [5, H]]...], (26) 

z — ' 777, 



where 



* = EE 



{W ^) 



AA (k)cr 



E 



AA 



Efc 



"A Cfe 



(27) 



Here (A;) stands for the electron or hole states secluded 
within a layer ±D around the Fermi level. E\\ = 
E\(D) — E\{D). The effective Hamiltonian with the 
charged states |A) = |lfocr), |36<r) frozen out can be ob- 
tained within first order in S. It has the following form, 



A (fc)cr 

/ , / , J kk> ** C fc<7 C fc'<T') 



AA'A kk'uu' 



where 



7 aa' 



^a'a 



JLv = (wrrw* 



i 



-Eaa — £fe E\> 



£fc' 



In the charge sector N = 2 the constraint J^ A A = 1 is 
valid. As is shown in pUlQ] , the effective Schrieffer- Wolff 
Hamiltonian of DQD describes not only the conventional 
indirect exchange between localized and itinerant spins. 
It also contains terms that intermix the singlet and triplet 
states of the quantum dot. This mixing is due to the tun- 
neling exchange with electrons in the metallic reservoir. 

As a result of integrating out the high-energy charge 
degrees of freedom, the effective spin Hamiltonian of the 
DQD acquires an 5*0(4) symmetry, which is the dynam- 
ical symmetry of a spin rotator. Before writing down the 
pertinent spin Hamiltonian a few words about a quantum 
spin rotator are in order. It is known (see, e.g., pq| ) that 
the symmetry of a standard quantum rotator is described 
by the operator of rotational angular momentum L and 
an additional vector operator M. These two operators 
generate the semi-simple algebra 04, they are orthogo- 
nal, L ■ M — 0, and the corresponding Casimir operator 
is L 2 + M 2 . The matrix elements of the operator M con- 
nect states with different values of the orbital momentum 
I — * I ± 1. The existence of this second operator reflects 
the "hidden angular symmetry" of the rotator. 

Similarly, the spin symmetry of the DQD is character- 
ized not only by the spin one vector S: one can introduce 
a second vector operator P orthogonal to S which deter- 
mines the matrix elements of transitions between the dif- 
ferent states of the rotation group 50(4). In the present 
case, the vector P = {P z , -P ± } determines the transitions 
between the singlet state and the different components of 



the spin triplet. It is convenient to express the spherical 
components of the vector operator P in terms of Hub- 
bard operators X AA (for brevity, AA' will be either S for 
singlet or /u = 1,0, 1 for the triplet magnetic quantum 
numbers), 



P+ = V2 (x ls - X s1 ) , P~ = V2 (x sl - X 



-is 



R = - (X as + X so ) 



(29) 



Within the same procedure, the spherical components of 
the spin one operator S are given by the following ex- 
pressions, 



5+ = V2 (A 10 + A 01 j , 

5- = V2 (a 01 + A 10 ) , S z = A 11 - A 11 . (30) 



The vector operators P and S obey the commutation 
relations of the usual 04 Lie algebra, 



(28) [Sj> &k] = tejkiSi, [Pj,Pk] = tejkiSi, [Pj,Sk] = icjkiPi 

(31) 

(here j, k, I are Cartesian indices). Besides, the following 
relations hold. 



S • P = 0, S 2 = 2(1 - X bb ), P 2 = 1 + 2X bb . (32) 

To wit, S (a spin 1) and P are two orthogonal vec- 
tor operators in spin space which generate the algebra 
04 in a representation specified by the Casimir operator 
g2 _|_ p2 _ ^ This justifies the qualification of DQD as a 
spin rotator. 

Returning back to the effective spin Hamiltonian, ( pq ) 
it now acquires a more symmetric form , 



H = H s + H T + H ST , 



(33) 



where 



H s = E S X SS + J S J2 X 



ss„ 



H T = EtJ2 X ^ + jTs ' s + —E *"""<" 

jjST = jST (p . s y (34) 



The local electron operators are defined as usual 
n a =c\c a =E c L c fc°-> s = 2 ~ 1/2 EE c ^ fcfe ' CT '' 

(35) 



kk' 



kk' aa' 



and f are the Pauli matrices. Moreover, the coupling 
constants are, 



jT= ( \Wl\ 



\W r 



\Ef — ei E r + Q — Ep 



(36) 



in case (b) and 



C. Two-spin representation 



J T = -\Wif 



1 



1 



v £f - ez Ei + Qi — ef 
in case (c). In both cases 



J s =a 2 ss J T , J ST =a ss J 7 



(37) 



(38) 



This completes the derivation of the spin rotator Hamil- 
tonian for a DQD hybridized with itinerant electrons. 



B. RG flow of coupling constants: Kondo 
temperature 

Due to the intermixing term H$t in the spin Hamilto- 
nian (f34|) both triplet and singlet states are involved in 
the formation of the low-energy spectrum of the DQD. 
Scaling equations for the coupling constants Jt, Jst can 
be derived by the poor man's scaling method of Ref . J24J . 
Neglecting the irrelevant potential scattering phase shift 
and using the above mentioned procedure of integrating 
out the high-energy states, a pair of scaling equations is 
obtained, 



dji 
dlnd 



[til 



{32? 



dn 



d\nd 



-2jiJ2. (39) 



(here_ ji = poJ T ,J2 = poJ ST ,d = p D). If S = 
Et{D) — Es(D) is the smallest energy scale, the energy 
spectrum of the DQD is quasi degenerate, and the sys- 
tem ( p9[ ) is reduced to a single equation for the effective 
integral j+ = j\ + j 2 , 



dj+ 
of In of 



= -ti+f 



(40) 



Then the RG flow diagram has an infinite fixed point, and 
the solution of eq. (Eoh gives the Kondo temperature 



T X0 = J Dexp(-l/j + ). 



(41) 



In the general case, the scaling behavior is more com- 
plicated. The flow diagram still has a fixed point at infin- 
ity, but the Kondo temperature turns out to be a sharp 
function of 5. In the case 5 < 0, |<5| ^> Tko considered in 
|7|,p|p0|] the scaling of J ST terminates at D ~ 6. Then 
one is left with the familiar physics of an under-screened 
S=l Kondo model |28| . The fixed point is still at infi- 
nite exchange coupling Jt, but the Kondo temperature 
becomes a function of 5. It is shown in Ref. S that a 
kind of universal law for Tk{5) exists also in this limit 



Tk/T K o = (T K0 /S)\ 
where 7 is a numerical constant. 



(42) 



It is known pq] that the algebra 04 can be represented 
as a direct sum of two 03 algebras. In our case this means 
that one can construct another pair of orthogonal oper- 
ators 



S+P S-P 

S l ^ n ) S 2 — 7: ■ 



(43) 



(see also ||). In the Hubbard representation the compo- 
nents of these spin vectors have a form 



S + = ^(X w + X 01 ± X 1S T X s1 ), 
' \/2 



S~ 2 = 4=(X 01 + X w ± X s1 T X 1S ), 



V2" 



S zh2 = -(X^+X^ T X V! >TX^). 
It is easy to check by direct substitution that 



(44) 



Sf = 3/4, X ss = - - (Sj • S 2 ), Y. X 



/,/! = i + (Si-s 2 ) 



(45) 

The Casimir operator can be introduced as 4S^ = 4S2. 

Then substituting (|43|) and (Q) in the Hamiltonian 
(|33| ) we rewrite it in the form 

H = J(S 1 ■ S 2 ) + J 1 (S 1 ■ s) + J 2 (S 2 ■ s) + 

J 3 (Si • S 2 ) 22 n<J + cons t ■ (46) 

Here 
J = Et — Eg = S, J1.2 = J ± J , J3 = -Jt — Js- 

(47) 

Thus, as was mentioned in Ref. |3|, the transformation 
( [43] ) maps the Hamiltonian (|3^) on an effective two-spin 
Kondo Hamiltonian plus an additional potential scatter- 
ing term. However, the physical meaning of these two 
spin operators differs from that in the conventional two- 
site Kondo model [g(|. They only span the two 03 sub- 
algebras of the semi-simple Lie algebra 04 (see next sec- 
tion for further discussion). 

This kind of effective Hamiltonian appears also in other 
situations where singlet and triplet states of a nanoob- 
ject are close in energy, e.g., in vertical quantum dots 
p,P| or in conventional dots at even occupation, pro- 
vided low-lying triplet excitons are taken into account 
H||. It was noticed in R|,0 that the interplay be- 
tween two energy scales, i.e. the interdot singlet-triplet 
gap 6 and the tunneling induced Kondo binding energy 
for triplet configuration A^ ~ Dexp(— 1/pqJt) results 
in essentially non-universal behavior of the Kondo tem- 
perature Tk (see preceding subsection). 



D. Magnetic field induced Kondo effect 

Yet another peculiar manifestation of " hidden symme- 
try" of the spin rotator is the possible occurrence of a 
magnetic field induced Kondo effect. Such possibility was 
discussed theoretically in Ref. j|] for the case of quantum 
dots formed in GaAs heterostructures and in Refs. |§|,||J7| 
for the case of vertical quantum dots where the exter- 
nal magnetic field influences the orbital part of spatially 
quantized wave functions and results in singlet-triplet 
level crossing. In DQD, similar effect arises if 6 > 
where the ground state of the DQD remains a singlet in 
spite of the tunneling induced renormalization. Here we 
re-derive the field induced Kondo effect in terms of spin 
rotator representation. 

In an external magnetic field, the energy levels in 
H T are split due to the Zeeman effect, Et — > -Etju = 
Et — [J.Sz- As was noticed in [9, the Zeeman split- 
ting Sz — 9^-bB of the excited triplet state compen- 
sates the energy gap 5 at a certain value of magnetic field 
B = Bq. In the vicinity of this point when 8 — Sz <C S 
only the levels Eti and E$ survive in the diagonal part 
E$X SS + E T XL X MM of the spin rotator Hamiltonian 
(p4J). Then the only renormalizable coupling parameter 
in the exchange Hamiltonian ( j34| ) is J ST . It is easily 
seen that the operators P + , P~ , (S z — X ss ) form an al- 
gebra 03 in the reduced spin space {S,T1}. In this sub- 
space the operators P + and P~ are reduced to \/2X ls 
and \/2X sl , respectively. The operators S + — v \/2X 01 
and S" — > \/2X 10 together with a new combination 
(A 00 -A n ) act in the subspace of excited states {TO, Tl} 
divided by the Zeeman energy from the low-energy dou- 
blet. These operators form a complementary algebra 03, 
and the direct sum of these algebras represent a realiza- 
tion of the SO(A) symmetry for a "spin rotator in an 
external magnetic field" when the rotational symmetry 
in spin space is broken. 

As a result, the effective spin Hamiltonian ( |34| ) in a 
subspace {S, Tl} reduces to 

H z = E z Ro + J ST (R • s) + H p . (48) 

Here Ez — E$ — Eti is the degenerate ground state 
energy level of DQD in magnetic field B = Bq, H p de- 
scribes irrelevant potential scattering, and the operators 
Ro, R are 



R = X 11 + X ss , R z = A 11 - X ss , 
R + = V2X 1S , R~ = V2X Sl . 



(49) 



The complementary vector T defined as 

T z = A 00 - A 11 , T+ = V2X°\ T~ = V2X 10 , (50) 

forms a second subgroup. This vector is quenched by the 
magnetic field. The spectrum of conduction electrons is 



also split due to Zeeman effect, but this splitting does 
not affect the Kondo singularity in the tunnel current: 
one simply may redefine the conduction electron energies 
and measure them from the corresponding Fermi levels 
for spin up and down electrons || . 

Applying the poor man's scaling procedure |2J] to the 
Hamiltonian (|48|), one comes to a scaling equation 



(1)2 



d\nD 



<hf 



(51) 



with a fixed point at j% — 00 and the Kondo temperature 
Tkz = Dexp(— I/J2), so that 



TjCQ 



= exp 



1 



1 



(52) 



Of course, the same kind of separation is possible for a 
degenerate pair of states Es, E T i, and the corresponding 
vectors R' and T' may be obtained from (|49] ) and (J5(j) by 
interchanging indices 1 and 1 (see Ref. |8j for a physical 
realization of this situation) . To summarize the descrip- 
tion of basic manifestations of spin rotator symmetry in 
DQD, we considered three limiting cases of a spin rota- 
tor representations depending on the physical situations: 
quasi degenerate state \5\ <C Tk when the resonance 
properties of a DQD are determined by the full SO (4) 
symmetry (|4l|), triplet ground state \S\ > Tk where the 
virtual excitations to singlet state render the Kondo tem- 
perature to be dependent on the initial singlet/triplet 
splitting ( |42j ) , and magnetic-field induced doublet ground 
state (^l|) where the Kondo resonance arises in spite of 
the loss of local rotational invariance. The hierarchy 
of Kondo temperatures is non-universal. The maximal 
value of Tk is given by Tko (§4]) , from which it falls with 
removing singlet/triplet degeneracy |7],|| . It then reaches 
the limiting value of Z)exp(— 1/j'i) at large S where the 
contribution of the high lying singlet state becomes neg- 
ligibly weak, and one returns back to the usual SU(2) 
symmetry of spin S=l described by the 03 algebra. 

All the above results could be obtained also in the rep- 
resentation (|44| ) . In this case the scaling equations should 
be derived for three coupling constants J\, J2, J3 of the 
spin Hamiltonian (Eq). This procedure is described in 
Ref. II . As expected, the results are equivalent since the 
scaling of J3 adds nothing to the singular behavior of the 
relevant parameters Jt and Jst- The problem becomes 
more complicated in case there are two sources and two 
drains pqlly] ■ Then, an additional index a should be in- 
troduced for the lead electrons, c^ a —> c a ka- The states 
with different a are intermixed due to interdot tunneling 
and one more operator S x P or Si x S2 should be intro- 
duced in the theory of spin rotator coupled to metallic 
leads il. 
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V. TWO-CENTER KONDO MODEL FOR DQD 

It is natural to expect that in the limit of vanishing 
interdot coupling V the tunneling through a doubly oc- 
cupied DQD is defined by the individual spins S = 1/2 of 
the left and right dot, and the existence or non-existence 
of resonance tunneling channel will be predetermined by 
the competition between the Kondo effect for the left and 
right dot and indirect exchange induced by the same tun- 
neling. In this limit the problem is reduced to a specific 
version of the two site Kondo model. The corresponding 
theory for symmetric two-site Kondo impurity in a metal 
was discussed, e.g., in Refs. [M. The question is, how 
this approach should be modified in the asymmetric cases 
(b) and (c). On the other hand, with increasing V, these 
two approaches should be matched, and it is instructive 
to compare the description of a DQD by two fictitious 
spins (EH) and by two real spins S; and S r . 

To study this problem we use the approach mentioned 
in Section 2, and consider the DQD as a two-site center 
with spins 1/2 in each site within the framework speci- 
fied by the Hamiltonian (m), m). In terms of Hubbard 
operators (10) for spin 1/2, this Hamiltonian is written 
as H = Hdo + H t + Hi r , where 



H do = J2 £ »A^f (A = 0, a, 2) 



i,A 

H t = ^\W fktr {Xf + Xf)c ka + h.c] . (53) 

In case (b), the states Ei , E\„ = ej, E riJ = e r and 
E r 2 = 2e r + Q are involved in the RG procedure, and 
the corresponding interdot tunneling Hamiltonian is rep- 
resented in the form 



H lr = V 



Z^ 



(XfXf + h.c.) 



(54) 



This tunneling is possible only in the singlet configura- 
tion of the DQD. We start with eliminating the polar 
states {10, r2} that arise due to interdot tunneling (54). 
This procedure, known as Harris-Lange canonical trans- 
formation |27[| , eliminates the interdot term Hi r , and in- 
stead, in second order in V, an interdot spin-Hamiltonian 
emerges, 



tllrs — J It 2^/1 r > 



(55) 



where Jy = T^ 2 /Ai. 

Like in the previous section, we should integrate out 
the high-energy charged states by using the Haldane's 
RG procedure p3| for the left dot alone, since the renor- 
malization of the deep level e r is negligible. However this 
procedure should now include the renormalization of Ji r 
due to reduction of the conduction band. The scaling 



equation for Ei is the same as (|24J) for a triplet state. We 
rewrite it in terms of the two-spin Hamiltonian as, 



dei _ Ti 
d\nD TT ' 



(56) 



(r; = Tt)- The same mapping procedure as in eq. ( JlSJ ) 
gives the following correction to the indirect exchange 
coupling constant 



Jir — Jir 



t%Ti\6D 
D ' 



(57) 



and its iteration results in the second scaling equation, 
dJi r 



d\nD 



-Pi-- 

7T 



(58) 



This procedure stops at D = D where the Schricffcr- 
Wolff limit for q is achieved. Integrating eq. (j^) from 
D to Do one comes to the following equation for the 
renormalized indirect exchange coupling 

J{D) = J{D )-(3lT l \n{D /D). 

Then, taking into account that 5 = Ji r by its origin, this 
equation can be rewritten in the form 



6 = S ~(3 2 l T l h-i(D /D). 



(59) 



This is the same result for the renormalized singlet-triplet 
excitation energy that we have found in the preceding sec- 
tion. In case (c) , a similar procedure starts with eliminat- 
ing the polar states generated by the interdot tunneling 
term (see fig. 3c), 



H lr = Y,\XfX^ + H.c]. 



(60) 



In this case the Harris-Lange procedure results in the 
indirect exchange Hamiltonian ( p5|) with coupling con- 
stant Ji r = V 2 J 'A2. Again, only the energy level ei is 
involved in the Haldane's RG procedure. The scaling 
equation (pq) contains on the right-hand side the factor 
/?2 instead of /3i, and its solution for the triplet/singlet 
level splitting gives 



S = 6o-^T l bi{D /D). 



(61) 



This is exactly the result obtained in Ref. JToj ] . 

Next, the Schrieffer-Wolff transformation eliminates 
the tunnel coupling. The operator S in (Eq) has the form 



*=£ 



Wi 



£k - £l 



(Xfb ka - H.c 



E 

ko 



Wi 



Zl+Ql- £k 



(Xf s b ka - h.c.) , 



(62) 



in case (b), and 
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* = £ 



Wi 



ka 



£k - £l 



(Xf°b. 



ka 



(63) 



in case (c). As usual, in second order, the tunneling term 
generates an indirect exchange between the leads and the 
dots. As a result, the total spin Hamiltonian acquires the 
form 

H s = J lr (Si • S r ) + J2 J '( S ' ' s ) + H '- ( 64 ) 

i— l,r 

Here J\ r = 5 is the renormalized singlet/triplet excita- 
tion energy ( p9[ ) or (|6l]) in cases (b) and (c) respectively. 
The components of the local spins Si are now defined as, 



5 



+ _ x n 



x 



it 



Ji 



1 



(X 



rr 



x 



iii 



(65) 



The Heisenberg-like interdot exchange (p4) arose in sec- 
ond order in V similarly to the effective AFM exchange in 
a half- filled Hubbard model [E7[. The indirect exchange 
coupling constants govern the interaction of the conduc- 
tion electrons and the local spins in the dots. They are 
given by, 



Ji = 



\Wtf 



Ep — £l 






\W r 



Ei+Q - E F ' 



(66) 



in case (b) and 



Ji = \Wif 



£f — £i £; + Qi — £f 



, Jr = 0, (67) 



in case (c). The last term H' in Eq. (J64|) includes irrel- 
evant potential scattering terms arising in second-order 
Schrieffer- Wolff transformation, and other invariants that 
appear in the Hausdorff expansion (26) i.e. mixed prod- 
ucts like Si • [Sj x s]. These terms arise due to interplay 
between the interdot exchange Hi r and the tunneling Ht- 
We do not consider here these small corrections to the 
main Kondo effect. 

Comparing (^4() with (flq), we see that the second or- 
der terms of the expansion reproduce the general struc- 
ture of a two-spin Hamiltonian. The interdot exchange 
has the same form for both representations, but there 
are significant differences in the values of the coupling 
constants between the leads and the local spins. In par- 
ticular, the tunnel coupling between the right dot and the 
leads is absent in this order in case (c), whereas in the 
two-spin representation ( |43]) the coupling constants are 
J\ = (2 + /?!) Ji and J^ =p^Ji- It should be noted, how- 
ever, that nonzero coupling between the leads and the 
right dot arises due to the interference between H\ r and 
S in a Schrieffer- Wolff representation, but its value dif- 
fers from J 2 - Thus, from the point of view of the general 
50(4) symmetry of the DQD, the two-site representation 
is simply one more representation of the 04 algebra as a 
direct sum of two 03 algebras. The only case when the 



representations Si, S2 and Si, S r coincide, is in the sym- 
metric DQD where the admixture of excitons is ignored 
and the parity is conserved for an isolated DQD. This 
symmetric DQD, of course, also obeys an 50(4) symme- 
try, but its ground state is a singlet. The only way to 
activate the "hidden symmetry" in this case is to switch 
on a strong magnetic field that compensates the exchange 
splitting. Then the effective spin Hamiltonian is given in 
(p8J) and a magnetic field-induced Kondo effect arises. 



VI. CONCLUDING REMARKS 

It is worth making several remarks about the ad- 
vantages of using alternative approaches to analyze the 
physics of the quantum spin rotator. We have presented 
three different ways of substituting spin 1/2 operators 
for the generic operators S and P of the 50(4) group. 
This substitution exposes numerous connections between 
the approach to the Kondo effect treating the double dot 
as a spin rotator and the conventional description as a 
two-site Kondo problem. The traditional theory of the 
two-site Kondo effect [gOj deals with a symmetric DQD, 
so it is formulated in terms of even-odd spin and charge 
states. The effects discussed in the present paper essen- 
tially arise only in asymmetric situations when 3\ 3> J% 
or J 1 3> J r . Besides, we treated conduction electrons 
in a single-channel approximation, whereas the even-odd 
state classification of conduction electrons in a two-site 
Kondo model pj| likes it to be a two-channel single- 
impurity model. Some generic properties of the two-spin 
Kondo effect are, nevertheless, similar in both limits. In 
particular, the competition between the on-site interac- 
tions Ji(Si ■ s) and the interdot exchange Jy(Sj ■ Sj) re- 
sults in the appearance of an unstable fixed point in the 
flow diagram dividing the Kondo singlet from an anti- 
ferromagnetic singlet ground states of the system. The 
conventional two-site Kondo impurity also can be classi- 
fied as a " spin rotator" , and the singlet-triplet (even-odd) 
mixing is an essential part of the Kondo physics in this 
case. 

Nano-objects whose symmetry is more complicated 
than localized spins, i.e. quantum rotors were discussed 
previously in a context of the theory of quantum phase 
transitions in 2D Heisenberg antiferromagnets, spin lad- 
ders and spin glasses (see, e.g. |30[). Quantum rotor was 
defined as a spin, whose rotation is constrained to move 
on a surface of M > 2 -dimensional sphere. An example 
of array of quantum rotors is a double layer system of 
antifcrromagnetically ordered quantum spins. If the in- 
terlayer coupling Kyi (Si,„ • S 2m ) in a site m is stronger 
than intersite coupling J mn (Si m • S ln ) in a given layer, 
the pair of spins Si m , S2 m form the quantum rotor S m 
with 03 algebra and Casimir operator 5 2 . The spin ro- 
tator is a natural generalization of this description: in 
case of M = 2 the excited triplet state in each site m 



12 



can be added to a manifold, and the ladder of spin rotors 
transforms into array of spin rotators. 

We leave more detailed discussion of electron trans- 
port through DQD for future communications, and dis- 
cuss briefly a limiting case of a biased DQD with J r = 0, 
which was realized experimentally [ojl^]. In case (b), 
this limit corresponds to an electrometer geometry . In 
this configuration the right dot is isolated from the leads, 
but, nevertheless, it can be used for driving the current 
through the left dot. In Ref. Q the driving was realized 
via electrostatic coupling between the dots, and charge 
transfer was allowed through monitoring the Coulomb 
resonance conditions. In the present case, the resonance 
in spin channel is allowed by modifying the energy of 
the charge transfer exciton. To measure this effect one 
should choose a E-window in a plan (V7,V') for sym- 
metric DQD. At zero difference VI — Vl no Kondo effect 
should be observed. Then, increasing this difference at 
given temperature T, one effectively changes the energy 



difference So (25) and rises the Kondo temperature Tk- 
When the regime T ~ Tk is achieved one finds oneself 
in a hatched region similar to that shown in Fig. 3a for 
a window {1,1}, and a zero bias anomaly should appear 
in conductance. 

In terms of the two-spin representation (|43j) , the struc- 
ture of RG scheme remains the same, and the only change 
is a disappearance of the second term in equation (aw for 
Jt- The changes in the real spin representation ( p5[) are 
more essential: J r — 0, and the Kondo tunneling occurs 
only through the left dot. However, the conventional the- 
ory of a single-site Kondo screening cannot be applied in 
this situation because the interdot tunneling term ~ J; r 
is still present in the Hamiltonian (64). If the renormal- 
ized coupling constant J; r remains positive, the left spin 
is dynamically screened and the right spin remains free. 
This model is a limiting case of under-screened spin-one 
solution. 

The description of the Kondo effect in terms of two 
fictitious spins R and T (49), ( pOJ ) is another example of 
separation of spin degrees of freedom into dynamically 
confined moment R and unscreened moment T (see also 
Rcf. || ), and this case was realized in experiments [QJ|] 
mentioned above. 

We have seen that the scaling trajectory for the cou- 
pling constant S = J; r is predetermined by the bare value 
of the singlet/triplet splitting So — 2/3iV, which, in turn, 
can be driven by the gate voltage [see eq. (jig)]. Thus, 
we see that the Kondo tunneling channel in the left dot 
can be opened by softening the charge transfer poten- 
tial that is governed by the right gate voltage Vg, and 
the DQD with isolated right dot works as a " charge-spin 
transformer" . 

The limit of zero J r in case (c) was considered within 
the spin-rotator approach in Ref. p0[ . In terms of a full 
50(4) description, the flow diagram is similar to that of 



a biased DQD (case b). If one would try to describe this 
asymmetric DQD in terms of screening of the individual 
spin S; , a problem would arise when taking into account 
charge fluctuations to the state \Exi) ([13]) at the first 
"Haldane" stage of the RG procedure, because this ex- 
citation is soft by assumption, uj ex j <C D. In this case 
the source of strong correlation is, in fact, the right dot, 
and the interdot tunneling is responsible for true spin- 
charge separation in the DQD. The description of Kondo 
tunneling in terms of the operators S and P is obviously 
preferable in this case. 

In conclusion, we considered here the spin excitation 
spectrum and the Kondo effect in DQD from the point of 
view of its generic symmetry, that is, an SO (4) symmetry 
of a quantum spin rotator. The properties of spin rotator 
differ in many cases from those of localized spin with the 
same 5 2 . In the case of triplet ground state (5 2 = 2) 
where cotunneling results in under-screened Kondo ef- 
fect, the presence of low- lying singlet excitation turns 
the Kondo temperature to be a non-universal quantity. 
If the ground state is a singlet (5 2 = 0), the Kondo effect 
is nevertheless possible if one projection of the low-lying 
triplet excitation is involved in electron tunneling. An al- 
ternative language for discussing the properties of DQD 
is provided by a two-site Kondo model approach. How- 
ever, in spite of the overall 50(4) symmetry of the prob- 
lem, the equivalence of these two approaches exists only 
in the case of conserved parity (symmetric DQD). When 
the asymmetric charge transfer exciton is admixed with 
a low-energy spin singlet, the two-site representation and 
the two-spin representation of the 50(4) group for a bi- 
ased DQD are not equivalent. 
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APPENDIX 

The wave- functions of symmetric and asymmetric 
DQD occupied by one two or three electrons are listed be- 
low, (see also p2fl). Besides, the tunnel matrix elements 
which connect the state from different charge sectors of 
the DQD are presented. The eigenvalues and eigenfunc- 
tions of an isolated neutral DQD with N = 2 can be 
found by direct diagonalization of the Hamiltonian (H), 
(H). In a neutral configuration {1, 1} the interdot capac- 
itive coupling is absent. Far from Coulomb resonances 
when the inequalities (flq) are valid, expansions ( J14| ) in 
symmetric case (a) and (|13|) in the asymmetric cases (b,c) 
give the following equations for the coefficients ay in first 
order of perturbation expansion in the tunnel coupling V. 
The processes taken into account in the mixing terms are 
shown by the dashed arrows in the upper panels of Fig. 
1. 



\3e,a) = j= (<$AA + dL4A) l°>> 



(a) Symmetric DQD: 



a\ 



V2/3. (A.l) 



(b) Biased DQD: 



o„ = 1 - &\ - P?,a s i = -oi. = V2/3[, 

a sr = —a rs = v2/3i, 

a ll = l-P' 1 2 ,a rr = l-Pl (A.2) 

Here /3[ = V/(Q + A). We assume that f3[ < (3 X and 
neglect the terms ~ /3[ in our calculations. 

(c) Asymmetric DQD: 



a ss = 1 -/?! - /3' 2 2 ,a a i 


= - 


(Jjgf — ^vs — V •"/■^2 1 




au =1—021 a rr = 1 - 


■0'i 



au 



V2/3 2 



(A.3) 



Here (3 2 = V/{Qr — A) -C /3a, and we neglect the corre- 
sponding contributions as well. 

To complete the enumeration of states involved in the 
tunneling Hamiltonian (Fll) one should define the charge 
states of DQD which arise in a process of electron tun- 
neling between the DQD and metallic leads. We are in- 
terested in one-electron tunneling, so the states with one 
and three electrons in DQD, N = 1,3, should be taken 
into account. 

(a) Symmetric DQD: 

One-electron states are even and odd combinations of 
electronic wave functions belonging to the left and right 
well. The same is valid for the three-electron states which 
in fact are the hole analogs of one-electron states. 

1 



(A.4) 



\le,a) = ^=[dl+dl)\0), 
l O) cr) = ±(dl-dl)\0), 



|3o,er) = 4= {dldl/ r] - dUMr) l°>- 



V2 

1 

(b,c) Asymmetric DQD: 

In this case the DQD is "polarized" both in negatively 
and positively charged states. The one-electron wave 
functions are the same in cases (b) and (c) 

\la,a) = (Vl - « 2 d\ a - a d\ a ) |0), 
|16, ( r)=(a4 + \/l-a 2 4 (T )|0} (A.5) 

(a = V/ IS). The corresponding energy levels are 

Eia = si + aV, E lb =e r - aV, (A.6) 

The three-electron wave functions are represented by ex- 
pressions 

|36, a) = (Vl - o? 44X T + a 4.«) |0), 



|3o, <t) = (-a d\jlA + V 7 !-" 2 4a4i d h) 1°) ( A - 7 ) 
in case (b), and 



|36, a) = (Vl - «' 2 dk« + a' d\J ri d^) |0), 

|3o, <r) = (-a' 4,« + V 7 ! - a' 2 4,4 A) 1°) ( A - 8 ) 

in case (c). Here a' — V/(Q r — Qi — A). The eigen-levels 
are given by the following equations 

E 3b = 2e r + Q - aV, E 3a = lei + Q + aV (A.9) 

in case (b) and 

E 3b = 2ei + Qi - a'V, E 3a = 2e r + Q r + a'V (A.10) 

in case (c). 

The tunneling matrix elements in the Hamiltonian H t 
(|ll| ) include states from different charge sectors N(i>i, v r ) 
of the dot Hamiltonian H c i (H). In the presence of 
an interdot coupling V, and at nonzero bias potential 
V 1 — V* > 0, the numbers vi, v r are non- integer, and the 
tunneling transparencies of the left and right dot are dif- 
ferent even if Wi = W r (case b) . In case (c) , the tunnel- 
ing barrier between the leads and the right dot is wider, 
and one can assume that W r < Wi, so that the symme- 
try is even stronger. Consideration of the asymmetric 
configurations in case (b), we note that the expansion 
coefficients in eq. (113) for the two electron states |A) are 
such that a ss 3> a s i, a sr (see eq. ( |A.2| )). The tunnel ma- 
trix elements which define the dominant contributions to 
the RG equations (EQ) arc 
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vv aucr /o 



w; 



ir 



f TO,lS 
T±,l± 
T±,3± 






v/2 



tO r , 



ir. 



5,1S 
qua 

S,3a 



1 

i 

71 



<7a ss Wl, 



dss^r 



(here u>; = yl — a 2 Wi). Similar equations can be de- 
rived in case (c) , where the wave- functions of the virtual 
charg ed states \lba) a nd \3b a) are given by eqs (A. 5), 
( |A~8| ). Now instead of ( |A.ll| ) one has 



W, 



W. 



TO,lS 
qua 

TO,3a 

qda 



T±,l± 



= -g W h W qu± 



w' t , W t 



qd± 



w. 



S,la 
qua 



V2 
-j=aa ss wi, W qda 



= wi, 
— aw r , 



1 

v/2 



where w', = yl 



/2 

a z wi. 
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